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Abstract 

Recent polarized proton experiments at Brookhaven National Laboratory are used as a basis 
for a model of the energy dependence of the analyzing power of proton-carbon elastic scattering. 
In addition to their practical value for polarimetry, the results of this analysis give constraints 
on the size of the Pomeron spin-flip coupling as well as information on the /a and u> spin-flip 
couplings. 
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During the past three years experiments have been carried out at Brookhaven National Labo- 
ratory, at the AGS with proton beam momentum of 21.7 GeV/c and at RHIC with beam mo- 
menta of 24 GeV/c and 100 GeV/c, to measure the single-spin asymmetry for elastic scattering 
of transversely polarized proton beams off (fixed) carbon targets. The principal goal of these ex- 
periments is to obtain a calibrated polarimeter based on elastic pC scattering in the CNI region, 
i.e. \t\ < 0.05 (GeV/c) 2 . In outline, the procedure to be used is to begin with an AGS beam of 
known polarization, and to measure the analyzing power A]y(t) near the RHIC injection energy. In 
the CNI region the analyzing power is expected to be a few percent which, coupled with the large 
rate, makes it well-suited to be a practical polarimeter. The ultimate goal is to have a sufficently 
precise calibration to make the desired 5% measurement of the polarization P. (As we will see 
the current experimental errors are too large to reach this level of absolute accuracy, though the 
relative errors are quite small.. This should improve soon as further experiments are carried out, 
but it may ultimately require the use of a polarized hydrogen gas jet target to reach the required 
level of 5% absolute). This AGS experiment is referred to as E950. See for a description of the 
experiment and the results. The next step is to measure the asymmetry after transfer of the beam 
into RHIC with a setup the same as or similar to the one in the AGS. The energy at this point is a 
little higher than where the calibration is made, but the analyzing power is not expected to change 
enough that significant additional uncertainty is introduced into the measurement of P. 

The difficult and as yet not solved problem has to do with the polarization at higher 
energy in RHIC. It would be nice to believe that the polarization is not reduced during the RHIC 
acceleration, but this is not really known, and it would be good to have some theoretical guidance 
as to what sort of energy dependence the CNI analyzing power is expected to have. The central 
goal of this paper is to present a plausible model, based on Regge phenomenology, to calculate 
A]y(t) at higher energy assuming only that a calibrated measurement can be made at injection. 
The model, as utilized here, has as spin-off some interesting theoretical results. In particular, the 
question of the spin-dependence ot the Pomeron coupling to the proton has been studied for many 
years with still no firm conclusion. (See [2] for an extensive discussion of this issue.) It is generally 
known to be small, less than 10-15%, but not better known than that. The analyzing power is very 
sensitive to this coupling, as we will see, and so even low energy data, like that from the AGS, can 
give us very good information about this coupling. 
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The model we will adopt uses three Reggeons, the Pomeron P and two Regge poles, f% and 
uj. Because carbon is an isosinglet none of the 1 = 1 Regge poles like the p or 0,2 can enter, but 
because they are known to make small contributions to the elastic, non-flip amplitudes we can 
use the fits for unpolarized pp and pp elastic scattering that have been carried out over an extensive 
energy range. |1] (To carry over our results to pp spin-flip will require important corrections 
because p and a2 spin-flip couplings are known to be strong [SJ E] - See Section 5, below.) 

We will begin by discussing the determination of A n (t) from the data obtained in E950. This 
analyzing power will then be applied to determine the polarization in RHIC from the asymmetry 
measured there at 24 GeV/c. Next we will describe how the 3-Reggeon model is extended from 
elastic non-flip scattering to spin-flip scattering and will examine what we can say about the model 
just based on the E950 results. We will then examine the RHIC data taken at 100 GeV/c and will 
see that, without knowing the polarization in RHIC at that energy, all the parameters of our model 
can be determined by measuring, in addition to the E950 results, the polarization-independent 
"shape" of the asymmetry distribution in the CNI region. This being done, we can calculate P at 
100 GeV/c without further data fitting. The model predicts the analyzing power at any energy, 
and we will examine this a little by checking on the variation between 21.7 GeV/c and 24 GeV/c 
and predicting the A]\[(t) at 250 GeV/c. The prediction for the 1 = part of the pp analyzing 
power through the RHIC colliding beam range will be given and some model dependence examined. 

We will attempt to include the p and a<i in our analysis in order to apply our results to pp 
scattering . The data from FermiLab E704 will be used to make tentative predictions for pp 
polar imetry at high energy. 

The effects we are looking at are rather small and the errors involved in the analysis are 
significant and generally correlated. Therefore we devote the last section of the paper to a careful 
analysis of the errors of the Reggeon spin-flip couplings determined earlier in the paper. The 
implications for the errors on the analyzing power as a function of energy are also worked out. 

2 

The formula we use for the pC analyzing power was derived by Kopeliovich and the author [2]. 
The approximations that go into the formula are discussed in that paper, and it is believed to be 
quite reliable over the small t-range involved in the CNI experiments. A particularly transparent 



2 



way to write it is 



= 1 - -Mr(s)] + -Im[r(s)]f(t), (1) 

with 

f{t) = ((1 + p 2 pC (t))(t/t c )(F£(t)/Fir(t)) - P P c(t) - 5 pC (t))/(l - Pp c(t)S pC (t)) (2) 

t(s) denotes the hadronic spin- flip parameter for the scattering. It will be defined precisely below. 
A]y(t, t) is the analyzing power in question and A]y(t,0) is "pure" CNI, the analyzing power in the 
absence of hadronic spin-flip. F^ m (t) is the electromagnetic form-factor and F^(t) is the hadronic 
form-factor for carbon; these are calculated in [S]. k = 1.79, t c = —d>^Za/a^ t and p p c(t) denotes 
the ratio of real to imaginary parts of the pC amplitude (It depends on t even if p for pp does not, 
as assumed in the derivation of Eq.l; it is also calculated in 8 .) 5 p c(t) denotes the Bethe phase 
(Hj for proton-carbon scattering; it is important only at the smallest values of \t\, smaller than for 
which data so far exists, but we carry it along anyway. We have neglected in writing Eq^ the 
dependence of the differential cross-section on r which is insignificant when r is small, in the range 
we anticipate. r(s) is defined by 



g 5 (s,t) = t(s)- — -g (s,t) (3) 
m 

go(s,t) and 55(5, t) denote, respectively, the spin independent and spin-flip pC elastic amplitudes. 
m denotes the proton mass. One of the main results of jS] is that r(s) is equal to the 1 = part 
of the corresponding spin-flip factor for pp scattering. It is, in general, complex and depends on 
s in an unknown way, but its t-dependence can be neglected over this small range of t. All the 
important t-dependence of Eq.l comes from the variation of the form factors and of p p c- The 
resulting dependence of f(t) is shown in Fig.l. 

It is important to note that this formula is given in terms of r(s) rather than 

f5°(s,t) =T(s)(i + p pC (s,t)) 

which is sometimes used. This has the advantage that, by the theorem of [Sj, t is independent of t 
while r\ has t-dependence inherited from p p c- 
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Figure 1: The t-dependent coefficient in the second term of Eq. 1. 



The E950 group fit their data in and determined values for rf . We could simply use 
their results, but because the propagation of errors to our later results is important, we have done a 
linear regression analysis of our own and calculated the error matrix for Re(r) and Jm(r). The data 
we fit is given in their Table 1 and we used the errors given there to determine the weights in the 
regression analysis. These errors include an approximately 12% error in the AGS beam polarization. 
(We have followed E950 in adding the errors linearly; perhaps this should be reexamined.) The 
results of this fits gives 



The very large error on the real part of t results from a modest error on the first term in 
Eq.l of about 25%,. The error in the imaginary part is also about 20%. These errors are too large 
for the job of doing a 5% absolute polarization measurement, and they will follow us through all 
of this work. 

There are additonal errors in these numbers because of unassessed errors in the parametriza- 
tion of Eq.l. In particular, the value of p p c used is calculated. It would be good to measure it. 

The same is true of the pC differential cross section which determines Fjj(t). The error in our 
calculation of these quantities is probably small, but it depends on, as input, the poorly known 
value of p for pp and pn at 21.7 GeV/c. This will need to be addressed when the other errors are 
reduced. (The values published in U agree within errors with Eq.Q.) 



t = (-0.214 ± 0.236) - (0.054 ± 0.015)< 



(4) 
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Figure 2: The error bands for the fit of E950 data to Eq.l 

The results of the fit are shown in Fig. 2 and Fig. 3, first showing the error bands on the 
regression analysis and second showing the Aj^{t) compared to the data. The error matrix for 



An(t) 




Figure 3: Fit of analyzing power using r = —0.214 — 0.054i to the data of E950 
Re(r) and Im(r) is given by 

2 f 0.0559 0.00332 \ 

° ' ' ~ \ 0.00332 0.000213 ) ^ ' 

and the corresponding 68.3% confidence level ellipse for the fit parameters {not the r's) is shown 
in FiglU 

It is important to bear in mind that if one attempts to fit the data using an incorrect 
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Figure 4: The la error ellipse for the fit to E950 data 



assumption regarding the polarization, call it P' , one will obtain a perfectly good fit with exactly 

the same \ 2 as t ne correct one; one will simply determine values of the fit parameters which scale 
as 



As a first application of these results, let us use the analyzing power just determined at 21.7 GeV/c 
to determine the polarization of the beam in RHIC after injection at 24 GeV/c. To do this we 
assume that A^(t) doesn't change significantly over this small energy range. We will return to this 
question after we develop a model for the energy dependence. We will use the preliminary data for 
the raw asymmetry e(t) which does not depend on any knowledge of the polarization of the beam, 
as presented at Prague in July 2002 JID a and at Spin 2002 [H]. The data for e(t) with errors e(t) 
are given in Table 1. The errors given assume that the systematic errors are equal to the statistical 
errors and they are combined in quadrature. (I thank D. Svirida for this information.) Evidently, 
these errors do not contain a contribution from the polarization because this is the raw asymmetry, 
not A^. (The values of An given in those references are calculated assuming P = 0.27 and so if 
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Im(r ) 



— Jm(r). 



(6) 
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one fits An instead of e(t) one will obtain values for r scaled as mentioned at the end of the last 
section. These turn out to be very different than the values in the E950 fit.) 

At this moment, the data is not yet published and must be considered in the context of 
conference presentations. In the not-too-distant future, new data will become available and this 
analysis may need to be revised. We choose to present this work based on the preliminary results 
for two reasons: (1) In order to illustrate the method we propose, it is necessary to apply it to real 
data, and (2) the results should be immediately useful to the ongoing program at RHIC in order 
to better understand the polarization at different energies. 



-t 


e(t) 


e(t) 


0.0117 


0.0064 


0.00060 


0.0138 


0.0059 


0.00043 


0.015 


0.0046 


0.00068 


0.0184 


0.0019 


0.00076 


0.0194 


0.0024 


0.00054 


0.0217 


0.0025 


0.00087 


0.0250 


0.0016 


0.00068 


0.0251 


0.0033 


0.00098 


0.0306 


0.00005 


0.00081 


0.0361 


0.00035 


0.00098 


0.0417 


-0.00046 


0.00116 


0.0474 


0.00089 


0.00138 



Table 1: -t, raw asymmetry e(t), and errors e(t) for RHIC 24 GeV/c 
Starting from the relation 

e(t) = PA N (t), 

we minimize the chi-square 

8=12 

X 2 (24) = £ (e(ti) - PA N (21.7, t t )f / e(t t ) 2 . (7) 

An (21.7, U) denotes the form in Eq.Q with r = —.0214 — .054i evaluated at ij. This is easily done 
with the result that 

P(24) = 0.39 ± 0.02. (8) 

The error here is just the statistical error of the fit. There is a larger error associated with the 
error in the determination of Apj(t). It can be determined using the formula for A^(t) and the 
error matrix, Eq.© to be about ±0.068 or about 17%. See Fig03 
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Figure 5: The uncertainty in the analyzing power predicted at 24 GeV/c 

Now using the E950 value of r and this polarization we plot together the expected asym- 
metry along with the RHIC 24 GeV/c data, Fig HO The fit is reasonable with x 2 = 1-5 per dof. 




P=0.39 



Figure 6: Raw asymmetry at 24 GeV/c as measured and as predicted by E950 analyzing power with 
P=0.39 
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Can we say anything about the energy dependence of Apf(t) based on the the E950 or RHIC 24 
data alone? The simplest plausible Regge model to describe the energy dependence has three 
Regge poles: the Pomeron, the /2 and the to. There have been at least two groups 011] which have 
produced fits to the extensive data for pp and pp elastic scattering over an enormous energy range 
which includes almost all of the RHIC energy range. The lowest energy, which includes our starting 
point, is marginal to the fits. Nevertheless we will use these models as a basis for our own model. In 
these fits the f2 and the u are always treated as simple Regge poles, but the more complex nature 
of the Pomeron is taken into account by allowing it to be alternatively a simple pole a little above 
J = 1 , or multiple poles at J = 1 in order to produce single or double logarithms in the asymptotic 
energy dependence of the total cross section. For this work we have chosen two representitive fits 
from P| 3 a simple pole at 1.0933 or a multiple pole at 1 designed to give a log-square growth. We 
will work through the first case first and then indicate how the second differs from it. 

We begin with the parametrization of pp elastic scattering given by Cudell et al [3], though 
one might do the same thing using other parametrizations such as that of Block et al [I] . Since it 
is known that the elastic, non-flip scattering is overwhelmingly 1 = exchange, even at 24 GeV/c 
[12.., we will assume the Regge couplings that they determine are for the 1 = families and so 
directly applicable to pC scattering. The form they assume for the forward amplitude is 



0o(s,O) =gp(s)+g f {s)+g u) {s) 



0) 



with 



9p(s) 




(10) 




ys~"(cot-(i-7/) -*) 




normalized that Im(go(s)) = cr tot (s). The values of the parameters given by them are 



e = 0.0933, 7] = 0.357, rf = 0.560, 



(11) 



X = 18.79, Y = 63.0, Y' = 36.2, 
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with X, Y, Y' in mb. 

Our model assumes that the spin-flip pp I = 



exchange amplitude g§(s,t) is given by 



firAt.l) - r(«)^g (s,t) (12) 
t {rpgp(s) +T f g f (s) +T u) g ul (s)}. 



m 



where r(s) depends on energy but not on t over the CNI range. It is in general neither real nor 
constant in s and is given by 

r(s) = {T P gp{s) +T f g f (s) +T u) g u (s)}/g (a,0) (13) 

where the Tj's are energy- independent, real constants. The phases of the amplitudes come only 
from the energy dependence as given in Eq.(|10jl. This is the key assumption from Regge theory 
which we need: as a result the real and imaginary parts of t(s) are given at each energy in terms 
of the three real constants Tp,Tf and t u . 

Since from E950 we can determine only two real parameters, it is plain that we can not fix 
this model completely. The best we can do is obtain a relation between the three spin flip couplings; 
for example, in FigQwe plot -ry and r w as constrained by the values of Eq.Q). 

The recent runs at RHIC have also provided data at 100 GeV/c. The polarization is not 
known there, so is it possible to obtain useful information from it? The answer is yes: if the 
asymmetry is described by the CNI formula, which is our fundamental assumption, then a fit to 
the raw asymmetry will determine (1 - |i?e[r(100)])P(100) and |Jm[r(100)]P(100). (In this and 
the following sections we will write the argument of the shape 5 and r as the lab momentum rather 
than the corresponding s. When we briefly discuss colliding beams, we will revert to s, but there 
should not be any confusion.) Thus the "shape" S of the distribution 

Im[T(p L )] 
k/2 - Re[T(pi)\ 

can be determined at any energy without knowing P. Now once we have the three quantities, 
Re[r (21. 7)], I m[r (21. 7)] and 5(100) we have three linear equations in Tp,Tp and r w and can deter- 
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Figure 7: Relations between the three Regge spinflip coupling based on E950 results alone 

mine them all. This then will give us t(pl) at pl = 100 GeV/c and thereby we can calculate the 
polarization P(100). We will now go through this process. 

Using data from the same sources as used at 24 GeV/c [TU1 ITT] . 



-t 


e(t) 


e(t) 


0.0117 


0.0036 





00055 


0.0138 


0.0029 





00047 


0.015 


0.0034 





00060 


0.0184 


0.0018 





00069 


0.0194 


0.0014 





00058 


0.0217 


0.0025 





00087 


0.0249 


-0.0004 





00072 


0.0251 


0.0009 





00085 


0.0306 


-0.0010 





00072 


0.0360 


0.0010 





00091 


0.0416 


0.0013 





00116 


0.0473 


-0.003 





00146 



Table 2: -t, raw asymmetry e{t), and errors e(t) for RHIC 100 GeV/c 



We then fit this as in Section 2 to the formula in Eq.Q with the right hand side multiplied 
by the unknown P(100). There is a small calculable energy dependence to f(t), and it is taken into 
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account in the fit. The result of the regression is 



2 

P(100)(l - ~Pe[r(100)]) = 0.263 (15) 
P(100)-Jm[T(100)] = -0.0137. 

K 

Combining these together we get for the shape of the distribution 

01 37 

5(100) = _^iJli = _o.052 (16) 
0.263 

By using Eq.@ the value of r(21.7) as determined in E950 in terms of the Regge spin-flip couplings 
via Ea, ()13|) and in the same way express S'(IOO) in term of the Regge spin-flip couplings via r(100) 
and Ea.l|13|). we have three equations to solve with the result 

Tp = -0.02 (17) 
T f = -0.43 
r w = 0.03 

There are significant errors in these determinations and we will return to them in the next section. 

These results allow us to calculate r(s) at any higher energy. (The model as it stands, is 
not really suitable for going to lower energy because lower lying Regge poles will rapidly beome 
important. Thanks to Boris Kopeliovich for emphasizing this limitation jS], |S].) Fig. |H] shows the 
pi dependence of the real and imaginary parts of r over the RHIC fixed target range. While the 
spin-flip couplings are evidently becoming smaller with increasing energy, there remains significant 
hadronic spin-flip at the top energy. 

We will use these results first for r(100) and thereby determine the polarization at 100 
GeV/c, P(100). From En.(fT3|) we find 

r(100) = -0.130 - 0.053i (18) 

With k/2 — Pe["r(100)] and/or im[r(100)], we can use the measured values of P(100)(k/2 — 
Pe[r(100)]) and/or P(100)7m[r(100)] to quickly determine P(100) = 0.23. Alternatively, as in 
Section 3, we can fit the 12 measured asymmetry values to An using Eq. ([18|) . Then minimizing 
the x 2 as hi Eq.fJTJ) we find 

P(100) = 0.23 ±0.02, (19) 
12 




Figure 8: Predicted s-dependence of r 

again the error is just the statistical error of this fit. In FigEJwe show the raw asymmetry measured 
at 100 GeV/c plotted with the prediction using Eq. (|18j) and this value of P. The agreement is 
reasonable, with \ 2 about 1.5/dof, about the same as the fit at 24 GeV/c. This is rather nice 
support for our approach. 

Given that P(24) was found to be 0.39, the value found for P(100) is surprisingly small. 
There seem to be three possiblilities to explain this: (1) There could be significant depolarization 
in the RHIC acceleration to 100 GeV/c. This is not expected, but this may be a signal of a problem 
in the acceleration. (2) The model used for energy dependence could be wrong. This seems very 
likely in its details, but if one compares the raw asymmetries at 24 and 100 GeV/c, one sees that 
a very large drop in the analyzing power-nearly a factor of 2 from 24 to 100 GeV/c is required if 
P remained constant. Of course, there might be a problem with the data at the smallest \t\, but it 
is not suggested by the errors assigned. 

In Section 3 we determined P(24) by assuming that t(24) is the same as r(21.7) determined 
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Figure 9: Raw asymmetry at 100 GeV/c predicted by model r(100) and polarization predicted to be 
0.23 

by E950. Now we have a prediction for t(24) and it is a little different: 

r(24) = -0.207 - 0.055?:. (20) 

We have checked its significance in determining P via Eq.(J7J) and find the best fit at P = 0.41±0.02, 
slightly different but certainly well within errors. (Remember that the result of the fits at 24 GeV/c 
were not used in the energy dependence calculations, so this is another check, although not a very 
strong one.) 

In the near future it is hoped to have a polarized proton run at RHIC at pi = 250 GeV/c. 
The predicted analyzing power is shown in Fig. 10. It is a little larger and has a slightly different 
shape from the 21.7 GeV determination. For clarity we also show in that figure the ratio of predicted 
to pure CNI analyzing power. This shows that the hadronic spin-flip must be taken into account 
in An even at the highest RHIC beam energy on a fixed target. 

As mentioned earlier, we also examined the behavior of a Cudell et al [3| fit with a log- 
squared asymptotic behaviour for the Pomeron. The fit we chose to look at has the Pomeron form 
replaced by 

gp(s) =i(A + Slog 2 s) + ttB logs, (21) 

s in units of GeV 2 . The Regge pole forms are unchanged, but the associated parameters are 
somewhat different: 

A = 25.29, B = 0.227, rj = 0.341, 7/ = 0.558, (22) 

Y = 52.6, Y' = 36.0, 
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Figure 10: Analyzing power predicted for 250 GeV/c pC scattering, and its ratio to the pure CNI 
analyzing power at 250 Ge V/ c 

with A, B, Y, Y' in mb. The spin-flip parameters associated with this form for the unpolarized 
amplitudes are found to be, in the same way as for the 3-pole model 



/ 

T P 



I 

Tp 



-0.02 
-0.49 
0.02. 



(23) 



These numbers are very similar to those found for the three pole model but not really directly 
comparable because the models are, in principle, quite different. They lead to very little difference 
in the model preditions. We can accentuate the difference by looking at the very high energy 
behavior where the difference between the power and the log's should be greatest, so in Fig llll we 
show the behaviour of Re[r] and 7m[r] over the full RHIC colliding beam energy range. (Note that 
here the x-axis is the s axis.) The difference remains tiny over the entire range. In either case the 
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1 = hadronic spin-flip should be negligible for 5% polarimetry for colliding beams above about 
70 GeV/c. 

-0 . 05 : 

Re(x) -o.i- 
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-0 . 02 
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-0 . 04 
-0 . 05 

100 500 5000 100000 S 

Figure 11: The 1=0 part of the pp r(s) plotted here against s through the RHIC colliding beam 
energy range. The solid line is the prediction of the 3-pole model, the dashed line the model with 
log-squared growth. 

5 

The last point leads us to make a brief observation. The r shown in Fig llll is only the 1 = part of 
the spin-flip, as predicted from the proton-carbon analysis. In order to use this for high-energy pp 
scattering we must add in the 1 = 1 spin-flip, in Regge terms the p and the 02, which are expected 
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to be quite large 0. If we fit the E704 data 012] we find 



t(E704) = 0.2 ± 0.2 + 0.024i ± 0.020i 

to be compared with our prediction of the 1 = part at pl = 200 GeV/c 

r (200) = -0.10 -0.0461 

Even with the very large error on the E704 result, these numbers are incompatible. The difference 
can be naturally explained by adding in large 1 = 1 spin-flip contributions which are absent in 
the pC scattering but are known to be important at lower energy in pp scattering Ej. If we 
assume that the p and have the same Regge behaviour as the lo and the /2, respectively, we can 
determine the C = — 1 and the C = +1 combined 1 = and 1 = 1 Regge flips to be, 



The relation of these numbers to the p and a2 couplings requires more study. However, just with 
these numbers and the value we have determined for rp we can calculate the values of the proton 
spin-flip expected at high RHIC colliding beam energy, and we find Tpp (70 2 ) = 0.04 + 0.02U and 
T pp (500 2 ) = -.01 + 0.006i so the hadronic spin- flip will have little effect on An in the CNI region 
at the highest energies. 



Finally, we would like to make some comments on how well the Regge parameters are determined 
in this model and using this data, especially on the question of the Pomeron spin-flip coupling. 
Starting from the error matrix for r(21.7) and error on S'(IOO), which requires the error matrix for 
r(100) we can easly propagate the errors to an error matrix of the 1 = Regge spin-flip factors. 
For completeness we give the error matrix for r(100) determined by our fit, and assuming the 
polarization P(100) = 0.23: 



T(C 



= -) = -2.10 



and 



T (C = +) = 0.69. 



6 




(24) 
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This leads to an error in the observed shape 5(100) = —0.052 ± .004. Using this plus Eq.(JSJ) leads 
to the corresponding la error matrix for the Reggeons 

/ T P Tf t u \ 

t p 0.0416 0.0559 0.0484 

T f 0.0559 0.0782 0.0612 

V t w 0.0484 0.0612 0.0618 / 

This may be visualized with the aid of three ellipses of projections of this ellipsoid onto the three 
planes defined by tp,tj and r w , in Fig. 




-0.2 -0.1 1 0.1 



X P 




Figure 12: Three projections of the error ellipsoid for the three 1=0 Regge couplings 



The errors are clearly very strongly correlated. Using these figures, one can respond to the 
question long asked |2j, is the Pomeron spin-flip coupling zero? Examination of these figures (or 
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the equation for the complete 3-d ellipsoid) shows that the point 



T P = 



(25) 



Tf = -0.44 



= 0.09 



is comfortably with all the ellipses. So the Pomeron can have a zero spin-flip coupling, but the 
/ cannot. Somewhat surprisingly, even at the highest RHIC energy and even in the absence of a 
Pomeron spin-flip, the pp spin-flip should not be expected to vanish although it is expected to be 
down at the few percent level. 

We can use the error matrix for the Regge couplings to calculate the errors in the predicted 
t(s) throughout the RHIC range. The errors vary slightly with energy but remain very close to 20% 
of either (tz/2 — i?e[r(s)]) or 7m[r(s)] between pi = 20GeV/c and pl = 250GeV/c. In particular, 
the error matrix implied by our model at pi = lOOGeV/c is 



We close by emphasizing that the numerical results obtained here are dependent on the 
preliminary results announced at Prague and Spin 2002 and they may change significantly. It will 
be very interesting to see how these results hold up when new, more precise experiments are carried 
out, both at the low end of the energy range here considered, and at the high end. It will also 
be important to extend the work of Section 5 to reach stronger conclusions regarding the 1=1 
couplings and to extend the Regge model downward in energy. The present level of accuracy on the 
experiments limits the strength of the conclusions we can draw, but it does appear that the CNI 
pC polarimeter should be a very precise relative polarimeter, with a modest absolute accuracy. 

I would like to thank Boris Kopeliovich, Gerry Bunce, Dima Svirida and Nigel Buttimore 
for very useful discussions of this physics. 




(26) 



larger than the errors of the fit Eq.lfH]). At pl 



250 GeV/c it becomes 




(27) 
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